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Discrete breathers, or intrinsic localized modes, are spatially localized, time-periodic, nonlinear 
excitations that can exist and propagate in systems of coupled dynamical units. Recently, some 
experiments show the sighting of a form of discrete breather that exist at the atomic scale in a 
magnetic solid. Other observations of breathers refer to systems such as Josephson-j unction arrays, 
photonic crystals and optical-switching waveguide arrays. All these observations underscore their 
in ■ importance in physical phenomena at all scales. The authors review some of their latest theoretical 

contributions in the field of classical and quantum breathers, with possible applications to these 
widely different physical systems and to many other such as DNA, proteins, quantum dots, quantum 
computing, etc. 
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^ i I. INTRODUCTION 

Oh 

Nonlinear physics of discrete systems has attained an enormous development in the previous years. A landmark is 
the prove of the existence of spatially localized, time-periodics oscillations in discrete nonlinear systems, by MacKay 
and Aubry |46|. They are now termed discrete breathers (DB), nonlinear localized excitations, or intrinsic localized 
modes, and ever since there are many experimental data and theoretical results published about these important new 
phenomenon in physics. They appear in lattices of oscillators interacting through nonlinear forces. The central class of 
system for which the concept of DB was developed is classical, autonomous spatially discrete Hamiltonian or reversible 
systems. A type of systems of this class is called Klein-Gordon lattices. They are Hamiltonian systems consisting 
of one degree of freedom anharmonic Hamiltonian oscillators coupled weakly in a lattice, that is, the frequency of 
oscillation of each oscillator varies non-trivially with the amplitude. Most studies on breathers have been done 
utilizing this type of models, with applications to biomolecules, such as DNA, solids, etc. Other classes of systems 
where breathers occur are: autonomous forced damped systems, and time-periodically forced systems. Considering 
applications to molecular crystals, quantum effects must be taken into account and the concept of quantum discrete 
breather should be developed. This review summarize some of the latest theoretical contributions in the field of 
localized nonlinear excitations carried out by the Nonlinear Physics Group of the Sevilla University. The structure of 
£h ■ this work has the following scheme: 
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• Section II: Breathers in Klein-Gordon lattices with competing short and long-range interactions 

— A: DNA model with dipole-dipole coupling 

— B: Numerical methods for obtaining breathers and determining their stability 

— C: Bifurcations 

• Section III: Moving breathers in homogeneous DNA chains 

— A: Moving breathers in DNA chains without curvature 

— B: Moving breathers in bent DNA chains 

— C: Parameters values of DNA 

• Section IV: Interaction of moving breathers with mass impurities 
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• Section V: Interaction of moving breathers with vacancies 

• Section VI: Small and large amplitude breathers in Fermi-Pasta-Ulam lattices 

• Section VII: Dark breathers in Klein-Gordon lattices 

• Section VIII: Polarons in biomolcculcs 

• Section IX: Quantum breathers 

— A: Quantum breathers in a translational invariant lattice 

— B: Trapping in lattices with broken translational symmetry 

• Section X: Thermal evolution of enzyme-created oscillating bubbles in DNA 



II. BREATHERS IN KLEIN-GORDON LATTICES WITH COMPETING SHORT AND LONG-RANGE 

INTERACTIONS 



Most of the studies about breathers in Klein-Gordon lattices refer to models which include only short-range 
interacting forces between the oscillators. This section reviews some of the results relative to breathers solutions in 
Klein-Gordon chains which include also long range interactions. 

A one-dimensional Klein-Gordon lattice is a chain of identical oscillators with an on-site potential V, each one of 
them coupled with its nearest neighbours through a parameter C <C 1. The Hamiltonian of a chain of N oscillators 
is: 



H = J2[ 2*» + V ^ + 2° {Un ~ Un - l)2 ) ' (1) 

where u n represents the coordinates of the oscillators referred to their equilibrium positions, and V(u n ) represents 
the on-site potential. 

In addition to the short-range interaction coupling, other types of forces can also exist. For example, if each 
oscillator has an intrinsic dipole moment, there exist also long-range forces due to the dipole-dipole coupling, and 
consequently the Hamiltonian must be modified. The existence of these two competing forces could alter substantially 
the properties of the breathers solutions. The study is presented in the next subsections, selecting a DNA model which 
has also into account the existence of dipole-dipole coupling. 



A. DNA model with dipole-dipole coupling 



In DNA there exists different kinds of interactions between the main atomic groups. One of them is the stacking 
interaction between neighbour bases along the DNA axis, these are short-range forces which stabilize the DNA 
structure and hold one base over the next one forming a stack of bases. There exist also long-range forces, due to the 
finite dipole moments of the hydrogen bonds within the nucleotides [3^, 0HJ . The dipole-dipole interactions between 
different base pairs become critical when the geometry of the double strand of DNA is taken into account, as the 
distance between base pairs and, therefore, the intensity of the coupling between them, depends on the shape of the 
molecule HEH^. 

The Peyrard-Bishop DNA model |5^| was initially introduced considering only the stacking interactions between 
base pairs. It can be modified by adding an energy term that takes into account the dipole-dipole forces. The 
Hamiltonian can be written as: 



H = ^2\ 2 mil n + v ( u n) + 2 C '( u »+i _ u «) 2 + 2 rr u «+P u « ■ (2) 

The term ^tom 2 represents the kinetic energy of the nucleotide of mass m at the nth site of the chain, and u n is the 
variable representing the transverse stretching of the hydrogen bond connecting the bases. The on-site potential is 
the Morse potential, i.e., D(e~ bUn — l) 2 , represents the interaction energy due to the hydrogen bonds within the base 
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pairs, where D is the well depth and represents the dissociation energy of a base pair, and 6 1 is related to the width 
of the well. The stacking energy is iC(u n _|_i — u n ) 2 , where C is the stacking coupling constant [5i| ■ The last term 
of the Hamiltonian is the long-range dipole-dipole interaction term, where J = q 2 /Aweod 3 is the dipole long-range 
interaction coupling constant |2l| . with q being the charge transfer due to the formation of the hydrogen bonds, and 
d the distance between base pairs, which is supposed to be constant. It is possible to consider D = 1/2, b = 1, m = 1 
without loss of generality. 

B. Numerical methods for obtaining breathers and determining their stability 

A static breather (SB) can be obtained by solving the full dynamical equations, which can be obtained applying 
Hamilton's equations to (|2J). This aim can be achieved using common methods based on the anticontinuous limit 
H, EH E3- US • The implementation of these methods basically consists in calculating the orbit of an isolated oscillator 
at fixed frequency wt>, and using this solution as a seed to solve the complete dynamical equations by means of 
a Newton-Raphson continuation method. A SB can be obtained for a given values of the dipole-dipole coupling 
parameter, J, and the stacking coupling parameter, C . This solution can be continued by varying C and maintaining 
J constant or vice versa. 

The linear stability of a breather can be studied by performing a Floquet analysis 10] . To this end, it is necessary 
to consider the evolution of a perturbation £ through the dynamical equations: 

in + V"(u n )Cn + C(2U - £ n+1 - &,_!) + V -£„ +p = 0. (3) 

p=£n 

The functions fi(0) = ({£„(0)}, {£„(0)}) must be integrated until t = T, giving fi(T) = ({£„(T)}, {L(T)})- There 
exists a matrix of dimension 2N called the monodromy, which is defined by fi(T) = A4Q(0). 

The spectrum of this matrix has the following property: if A is an eigenvalue, then A*, 1/A and 1/A* are also 
eigenvalues. Therefore, stability implies that all the eigenvalues have modulus unity 20]. 

The procedure of continuation of SBs can be done as long as an integer multiple of their frequency does not resonate 
with any of the frequencies of the linear modes (or phonons). This resonance provides an analytical expression for 
the upper boundary delimiting the existence of SBs: C — (uib/uj ) 2 — 1/4 + 3£(3) J/8. 

C. Bifurcations 

The study of the bifurcation loci of SBs is interesting for three different reasons: a) It gives the range of existence 
of the SBs; b) It gives the regions in the parameter space where the breathers are stable (and, therefore, physically 
observable in real systems) or unstable; c) Finally, but of paramount importance, it gives the values of the parameters 
where the breathers are movable. 

Two types of bifurcations appear as C and J are varying: 1) Stability bifurcations: they occur when a breather 
changes its stability. The existence of this type of bifurcation is a necessary condition for the existence of movable 
breathers (see Section llllfl . 2) Breather extinctions: they occur when a breather is not continuable any longer. 

When only stacking interaction is considered, only stability bifurcations occur and moving breathers are possible. 
However, with only dipole-dipole interactions there are no stability bifurcations and, therefore, moving breathers do 
not exist. 

When both interactions exist simultaneously, the analysis of 1-site (site-centered) and 2-site (bond-centered) 
breathers, permits to establish the following staments: A) 1-site breathers can be moved only for J < J c \. B) 
2-site breathers are movable for J < J c2 being J c2 > J c %. They are also movable for J S [J C 2, Jo)- Therefore, J 
establishes the maximum value of J (and also the minimum value of C) for which breathers can be movable, i.e., there 
are no moving breathers for C < J and for J > J a . This is the reason why J a is called the mobility limit pl| . 

The bifurcation loci for the stability bifurcations are shown in Fig. ^a). Also, in Fig. ^b) the dependence of the 
mobility limit with respect to the breather frequency is shown. This curve fits very well to the relation J D = Ato^, 
where A=0.1921 ± 0.0002 and r=2. 377 ± 0.005. 

III. MOVING BREATHERS IN HOMOGENEOUS DNA CHAINS 

This section is dedicated to the study of travelling breathers, or moving breathers, in DNA chains described by 
the Peyrard-Bishop modified model. The first subsection presents this study for straight line DNA chains, or chains 
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(a) (b) 




FIG. 1: (a) Bifurcation loci for the first and the second exchange of stability in 1-site and 2-site breathers. The values 
corresponding to lj^ = 0.8 are: J c i = 0.076, J C 2 = 0.085, J = 0.11263. (b) Dependence of J with respect to the breather 
frequency. 

without curvature. The next subsection considers DNA chains which have been bent and adapted to a parabola. 

A. Moving breathers in DNA chains without curvature 

Contrary to static breathers, moving breathers are not exact solutions of the dynamical equations of the model 
and they have a finite life due to the effects of phonon radiation. A moving breather is obtained by 'kicking' a static 
breather, i.e., perturbing its velocity components. However, not all breathers can be moved, and not every kick 
can move a breather. Therefore, it is important to know the conditions that a breather must fulfill in order to be 
movable and what the characteristics of the perturbations are. As indicated in 0, UtI , the following steps must 
be performed in order to obtain a moving breather: 1) To look for the existence of the two complementary stability 
bifurcations for the 1-site breather and the 2-site breather. Their bifurcation loci must have a region where they are 
fairly close. The static breather that can be moved should be obtained for values of the parameters close enough to 
these bifurcation loci. However, as shown in |2J it is only necessary that the static breather is "quasi-stable", i.e., 
not very stable or not very unstable; 2) To perturb the breather with the velocity components of the marginal mode 
at the neighbouring bifurcation. Time evolution of a moving breather is shown in Fig. 

An useful concept used for describing the breather dynamics is its effective mass [ll|, [lTj ■ It is a measure of the 
breather inertia to external forces. The breather velocity must be perturbed in a direction colinear to the marginal 
mode in order to obtain moving breathers. This mode is the responsible for the stability bifurcation. If the normalized 
marginal mode is V, the perturbation added to the breather velocities, which are zero at t = 0, is given by AV, and A 
is the magnitude of the perturbation. Thus, the kinetic energy added to the breather by the initial kick is 1/2 A 2 . It is 
found that the resulting translational velocity of the breather, v, is proportional to A. Then, the concept of effective 
mass can be defined through the relation m*v 2 /2 — A 2 /2. Therefore, m* = (X/v) 2 . Consequently, moving breathers 
can be considered as a quasi-particle with mass m*. The effective mass is a quantitative measure of the breather 
mobility. Larger mass indicates smaller mobility. 

Numerical studies show that, for high values of the dipole-dipole coupling, the predicted range of existence of moving 
breathers decreases and, in some cases, it is impossible for the breather to be moved even in the neighbourhood of 
the stability bifurcations. Another important result is that the dipole-dipole interaction affects to the mobility of the 
breathers. The study of the dependence of the breather effective mass, shows that the mobility decreases when the 
intensity of the dipole-dipole coupling increases (see Fig. |3J). 
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FIG. 2: Evolution of a moving breather for C = 0.1244, J = 0.06 and an initial 'kick' A = 0.10. The breather involves basically 
three sites that oscillate in phase. 



B. Moving breathers in bent DNA chains 

Suppose a DNA chain lying in a plane which has been bent and adapted to a parabola of curvature k (that is, 
the location of the n-th base pair is determined by the equation y n — kx^/2) |5j- All the dipole moments are 
perpendicular to the plane and parallel among them. Thus, the modified Hamiltonian of Eq. J5J is: 



H = ^2\ o mil l + v ( u n) + ^C(u n+1 - u n f + - |^ _ z l3 u n+p u n , (4) 

n=l \ Z 1 1 p^n |7 ™ ^ / 

where f n = (x n ,y n ), and it is assumed that the chain is inextensible, so that the distance between neighbouring sites 
remains constant: \f n — r n +i\ = d. 

Fig. ^a) illustrates the evolution of the energy centre [2lJ of a moving breather in a bent chain. If the added 
kinetic energy, E = A 2 /2, is smaller than a critical value E C) the breather rebounds, but, if E > E c , the breather 
passes through the bending point. Fig. E{b) shows that the critical energy increases monotonically with the curvature. 

Chain bending acts as a hindrance for the movement of breathers. This hindrance resembles to the experimented 
by a particle moving in a potential barrier. In this case, breather can be consider as a quasi-particle and this barrier 
can be calculated by finding the points where breathers rebound, i.e., the turning points for different values of E. 
Furthermore, if the breather has a constant effective mass m*, this potential barrier can be obtained using the 
expression: 



E b (x) = -\ 2 [l-(v(x)/v )% (5) 

where v(x) is the translational velocity and v Q is its value at t = 0. Fig. [S]show that there exists a good agreement 
between the barriers calculated using both methods for a given value of k. The barrier calculated by the second 
method exhibits an irregular shape, whose origin lies in the non-uniform behaviour of the translational velocity due 
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FIG. 3: Breather mass versus stacking coupling parameter C for dipole-dipole coupling parameter J = (left), J = 0.03 
(center) and J = 0.05 (right). Note that the maximum value of the mass occurs at the first bifurcation, and that the breather 
mass for J = 0.05 is always greater than the mass for J = 0. Note also the different scales at the Y-axis. 
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FIG. 4: (a)Evolution of the breather energy centre (Xe) for two different initial velocities. The curvature is k = 2. (b) Critical 
kinetic energy E c versus curvature n. 
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FIG. 5: Potential barrier calculated finding the turning points (solid line) and using Eq. © for E = 0.0162 (dashed lined) and 
E = 0.0200 (dotted line). The curvature is k — 2 and the critical energy is E c ~ 0.0126. The zero value of the energy centre 
(Xe) represents to the bending point. 



to the discreteness of the system [l9j. This result confirms that, in this case, a moving breather behaves as a particle 
of constant effective mass m*. 

It is worth remarking that a similar study has been performed in bent alpha-helix proteins Q ■ 



C. Parameters values of DNA 



The results presented are valid for all chains with both short range and long range interactions given by equations 
(J2J) and ffl. When they model a DNA chain, the following parameter values should be used [25|]: D — 0.04 eV, 
b = 4.45 A" 1 , d = 3.4 A, J = 0.0031 eV/A 2 , C = (0.01, 10) eV/A 2 , m = 300 amu and the time unit is 0.88 ps. 
For this set of values, the minimal kinetic energy for a moving breather to cross the bending point is small, and, in 
consequence, it seems plausible that most MBs will have enough energy to pass through bending points, although the 
ones with low energy will be reflected. 



IV. INTERACTION OF MOVING BREATHERS WITH MASS IMPURITIES 



Sometimes a chain can have a single point inhomogeneity, that is, the chain is homogeneous except at a single site 
where, for example, an impurity or a different type of particle is located. Moving breathers can be generated far apart 
the impurity and be launched towards it. This section presents some results relative to the effects of the collisions of 
moving breathers with an impurity in a Klein-Gordon chain. 
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FIG. 6: (a) Frequencies of the linear modes versus the parameter a. At a = a r es and a — a c , two different bifurcations occur, 
being the first one due to the resonance between the impurity mode and the breather, (b) Different regimes in the interaction 
of a moving breather with an impurity. 



Consider a Klein-Gordon chain with nearest neighbours attractive interactions, with Hamiltonian given by: 



H = J2 Q«n + KK) + \c{u n - u n ^) 2 ^j , (6) 

where V n (u n ) = D n (e~ Un — l) 2 is a substrate potential at the n-th site. The inhomogeneity is introduced assuming 
a different well depth at a single site, i.e., D n = D a (\ + aS n fi), where a G [— l,oo) is a parameter which tunes the 
magnitude of the inhomogeneity. The particle located at n — is an impurity. 

The dynamical equations can be linearized if the amplitudes of the oscillations are small enough. These equations 
have N — 1 non-localized solutions (linear extended modes) and one localized solution, (linear impurity mode). Their 
frequencies, loe and wl, are given, respectively, by: 



u>(q,a) = ^2 + 4Csin 2 ^, u> 2 L = uj 2 q + 2C + sign(a)y/oJ^ + 4C - 2) (7) 

where q £ (0, 7r] if a < and q £ [0, n) if a > 0. Fig. EJa) shows the dependence on a. 

Numerical simulations show that, depending of the critical values of the parameter a, there exist four different 
regimes in the moving breather-impurity interaction |2 ll |: 1) Barrier. The impurity acts as a potential barrier. It 
occurs either with a > or a £ (— 1, a{) with ct\ < 0. If a > 0, the breather can pass through the impurity provided 
the translational velocity is high enough 23] . Excitation. The impurity is excited and the breather is reflected. This 
case appears for a £ (011,012), an example can be seen in Fig. E[a). 3) Trapping. The breather is trapped by the 
impurity. It occurs in the interval a € (02,0:3). When the moving breather is close to the impurity, it becomes 
trapped while its center oscillates between the neighbouring sites, as Fig. Efb) shows. 4) Well. The impurity acts 
as a potential well. It occurs for a € (0*3, 0) and consists of an acceleration of the breather as it approaches to the 
impurity, and a deceleration after the breather has passed through the impurity. 

It is observed that the breather bifurcates with the zero solution at a = a res . That is, for a smaller than this value, 
no impurity breather exists. At a = a res , the frequency of the impurity mode coincides with the moving breather 
frequency, i.e., in 0, — ^b- 

The scenario for the trapped breathers when a < is the following: the impurity mode has q = 0, and also all 
the particles of the impurity breather vibrate in phase; this vibration pattern indicates that the impurity breather 
bifurcates from the impurity mode and it will be the only localized mode that exists when the impurity is excited for 
a > a res . Thus, when the moving breather reaches the impurity, it can excite the impurity mode. For a < a res , the 
moving breather is always reflected. In addition, the impurity breather does not exist. Therefore, there might be a 
connection between both facts, i.e., the existence of the impurity breather seems to be a necessary condition in order 




FIG. 7: (a) Interaction of a breather with an impurity for a — —0.52, which corresponds to the impurity excitation case, (b) 
Evolution of the moving breather for a = —0.3, which corresponds to the trapping case. The moving breather becomes trapped 
by the impurity; afterwards, the breather emits phonon radiation and its energy centre oscillates between the sites adjacent to 
the impurity. 



to obtain a trapped breather. If a > 0, the impurity mode has q = ir but the impurity breather's sites vibrate again 
in phase, that is, the impurity breather does not bifurcate from the impurity mode. There are two different localized 
excitations: the tails of the (linear) impurity mode and the impurity breather. Thus, if the moving breather reaches 
the impurity site, it will excite these localized excitations. Therefore, we conjecture that the existence of both linear 
localized entities at the same time may be the reason why the impurity is unable to trap the breather when a > 0. 

All these facts allows the formulation of the following trapping hypothesis: The existence of an impurity breather 
for a given value of a is a necessary condition for the existence of trapped breathers. However, if there exists an 
impurity mode with a vibration pattern different from the impurity breather one 's, the trapped breather does not to 
exist. 



V. INTERACTION OF MOVING BREATHERS WITH VACANCIES 



Another type of local inhomogeneity that can exists in a homogeneous chain is a vacancy, that is, the absence of 
a particle at a single site. A moving breather can be launched towards the vacancy and the effects of the collision 
should be analyzed. 

Consider a Hamiltonian Frenkel-Kontorova model with an anharmonic interaction potential |l6j : 



H = E 2*» + 4^2 t 1 ~ cos(2ira)] + C ^[exp(-6(x„ +1 - x n - 1)) - l] 2 , (8) 
n ~ 

where {x n } are the absolute coordinates of the particles. The choice of a periodic potential makes easy to represent 
a vacancy. Thus, if the vacancy is located at site n v (see Fig. [SJl, the displacements of the particles with respect to 
their equilibrium positions are: 

u n = x n — nh n < n v 

(9) 

u n = x n — (n + 1)L n > n v . 



It is possible to generate a moving breather travelling towards the vacancy at the site n v , and simulate the collision 
with this type of point inhomogeneity. 

The initial perturbation, {V„} has been chosen as V = A(. . . ,0,— l/v2,0, l/y/2,0, . . .), where the nonzero values 
correspond to the neighboring sites of the initial center of the breather. 
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FIG. 8: Scheme of the Frenkel-Kontorova model with sine-Gordon on-site potential, with full circles representing to the 
particles. 



When a moving breather reaches the site occupied by the particle adjacent to the vacancy, i.e., the location n v — 1, 
the particle can jump to the vacancy site or remain at rest. If the former takes place, the vacancy moves backwards. 
However, if the interaction potential is wide enough, the particle at the n v + 1 site, can feel the effect of the moving 
breather at the n v — 1 site and it can also move backwards the vacancy site. In the last case, the vacancy moves 
forwards. 

Numerical simulations show that the occurrence of the three different cases depends highly on the relative phase of 
the incoming breather and the particles adjacent to the vacancy |22|. However, some conclusions can be extracted: 
1) the incident breather always losses energy; 2) the breather can be reflected, trapped (with emission of energy) or 
transmitted by the vacancy, in analogy to the interaction moving breather-mass defect |24| : 3) the transmission of 
the breather (i.e. the breather passes through the vacancy) can only take place if the vacancy moves backwards, i.e. 
the particle to the left jumps one site in the direction of the breather. An explanation of this fact is that the particles 
to the right of the vacancy, in order to support a moving breather, need a strong interaction which cannot be provided 
by the interaction across a vacancy site, because the distance correspond to the soft part of the Morse potential. Fig. 
illustrates these phenomena. 

The moving breather-vacancy interaction is very sensitive to the velocity of the breather and the distance between 
the breather and the vacancy. Consequently, a systematic study of the state of the moving breather and the vacancy 
after the interaction cannot be performed. Therefore, we have performed a great number of simulations each one 
consisting in launching a single breather towards the vacancy site. In particular, we have chosen 1000 breathers 
following a Gaussian distribution of the perturbation parameter A with mean value 0.13 and variance 0.03 for different 
values of the parameters b and C. Fig. HOf a) shows the probabilities that the vacancy remains at its original site, 
or that it jumps backwards or forwards, for C = 0.5 and C — 0.4. From this study, an important consequence can 
be extracted: concerning to the parameter 6, there exist two different regions of values, separated by a critical value 
bo(C), such that for b > bo(C), the probability that the vacancy moves forwards is almost zero, whereas for b < 6rj(C), 
this probability is significant. 

The non-existence of forwards vacancy migration can be explained through a bifurcation. These bifurcations are 
related to the disappearance of the entities we call vacancy breathers. They are defined as breathers centered at the 
site neighboring to the vacancy, e.g. the n v — 1 or n v + 1 sites. It can be observed in Fig. HOf b) that, for b below the 
bifurcation value, vacancy breathers do not exist. 
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FIG. 9: Energy density plots of the interaction moving breather-vacancy. The vacancy is located at n v = 0. (a) the vacancy 
moves backwards; (b) the breather is transmitted by the vacancy; (c) the breather is reflected and the vacancy remains at rest; 
(d) the vacancy moves forwards. 



VI. SMALL AND LARGE AMPLITUDE BREATHERS IN FERMI- PASTA-ULAM LATTICES 

This section is dedicated to the study of breathers in Fermi-Pasta-Ulam (FPU) lattices, they can be considered as 
a subclass of the more generic Klein-Gordon lattices, in which the on-site potential term has been set to zero. An 
FPU lattice can be simply defined as a one-dimensional chain of particles subjected to nonlinear interactions between 
them. Therefore its Hamiltonian is given by 



-x 2 n + W(x n - z n _i) 



(10) 
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FIG. 10: (a) Probability that the vacancy remains at its site (squares), moves backwards (circles) or moves forwards (triangles), 
for a Gaussian distribution of A as a function of the inverse potential width b. (b) Amplitude maxima of a vacancy breather 
versus b. The vacancy breather disappears at the bifurcation point. This is related to the vanishing of the forward movement 
probability. 



where x n describes the displacement of a particle from its equilibrium position, and W(x) is a nonlinear interaction 
potential satisfying W(0) — W'(0) — 0, W"(0) > 0. A characteristic of this system is that it docs not possess the so 
called "uncoupled limit", for which there exist trivial breathers. Although some early results on DBs were obtained in 
FPU chains in the late eighties by Sievers and Takeno j5(| , a complete and rigorous proof of their existence in these 
type of lattices was provided only some years ago. In 2001, using a variational method, Aubry et al. 01 proved 
the existence of breathers with frequency above the phonon spectrum if the non linear interaction potential, W. is 
a strictly convex polynomial of degree 4. In the same year and based on a center manifold technique, James |40j 
proved the existence (respectively, nonexistence) of small amplitude breathers (SAB) with frequency slightly above 
the phonon band if W satisfies (respectively, violates) the following local hardening condition: 

B = iw w (0)-(W (s >(0)) a >0. (11) 

In spite of the implicit character of both proofs, the center manifold method provides additional information on the 
breathers amplitudes. In fact the method shows that for the interaction force y n — W'(u n ), all small amplitude 
time-periodic solutions with frequency Wb of the FPU system (including breathers) have the form 

y n (t) = b n cos(w b t) +h.o.t, (12) 

being the amplitudes b n determined by the two-dimensional map 

b n+1 + 2b n + 6 n _i = -fi b n + Bbl + h.o.t, (13) 

provided <x>b lies near the upper edge, u^, of the phonon band, i.e. provided /i = u> 2 — to 2 . <C 1. For B > fixed 
and ji > small enou gh, DBs correspond to homoclinic solutions to of the recurrence relation (|13fl satisfying 



lim b n = 0. In Refs [4(j and 41] an existence proof is given for two homoclinic solutions to 0, denoted as ±b n , 



ib n , having different symmetries &L„_i = — b n (bond-centered mode) and b 2 _ n = b n (site-centered mode). 

Since </i< 1 Eq. i|13|) can be approximated at leading order by an integrable differential equation [60| and it 
is easy to find the following approximations to the exact solutions y n and y n : 

„ i ,, i - ,_!,». zlL C0SUJht (14) 

UnK ' ~ 1 M/ B cosh( (|n + 1/2| - 1/2) «JJi) ' 1 ' 




v 2 (t) ~ ( 1)V 2M COSLJht (15) 
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FIG. 11: Comparison between a small amplitude breather (squares) for B > and a large amplitude breather (full circles) 
for B < with the same frequency Wb = 2.01. This last breather is superposed to an uniformly static state given by 
linin^ioo Un ~ 0.03. The dashed line represents the centre manifold approximation l|15t for small amplitude breathers. 
Parameters: SAB: K 3 = -0.3, Ka = 1 (B = 2.64); LAB: K 3 = -1, Ka = 1 (B = -1). 



These expressions show that the maximum amplitude of SAB is 




while their width is 0(^ _1 / 2 ) and diverges as Wb — > o>+. 

It is important to check the range of validity of the center manifold approximations (|14fl . I|15|) . A detailed numerical 
study of DBs in an FPU chain has been done [6(| with the anharmonic potential 

W(u) = £ + f « 3 + f - 4 (17) 

and fixing K4 — 1. With this choice _B = 3 — 4if|, and thus the parameter B is positive if \K$\ < if| = y/3/2 ~ 0.86. 
In the case = the potential is even with _B = 3 > 0. For numerical convenience it is better to use as dynamical 
variables the difference displacements u n — x n — x n -\- With these variables the equations of motions become 

u n + 2W'(u n )-[W'(u n+l ) + W'(u n - 1 )}=0 , neZ. (18) 

Note that u n is one-to-one related to the forces y n at small amplitudes, because W' is locally invertible since W"(0) ^ 0. 
Indeed, expressions (|14l) . I|15|) also approximates u n for small [i since u n = y n + 0(y„). Moreover ussing periodic 
boundary conditions w n +2p(i) = u n (t) the maximum frequency of the linear phonons, u x , is exactly 2 as in the infinite 
lattice. 

As can be appreciated in Fig. 1111 there exist an excellent agreement between the approximation (|15fl (dashed 
line) and an exact site-centered mode obtained numerically (squares) for c^b = 2.01 (/1 w 0.04) and K3 = —0.3. 
Continuing this SAB as a->b goes away from the phonon band, according to (| 1 fc>l> . the maximum amplitude of the 
relative displacements is approximately a linear function of fi 1 / 2 up to /1 ps 2, relatively far from the parameter 
region in which the center manifold approximation is valid in principle. In fact, the results show that, surprisingly, 
expressions (|14|l and l|15|) fit very well the profile of the relative displacements u n even for breathers with moderate 
amplitudes. 

The linear stability of the SAB can be studied through a Floquet analysis [lj} ■ I n the symmetric potential case the 
result is well known ^^|: the Page mode (site-centred mode in the u n variables) is stable while the Sievers-Takeno 
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mode (bond-centred mode in the u n variables) has a harmonic instability (a pair of real eigenvalues a, c -1 close to 
1) that increases with the breather frequency. When a cubic term in the potential W is introduced (K 3 ^ 0), the 
situation is more complex. Non-even potentials induce oscillatory instabilities in both modes, which increase with 
1 . Although most of the oscillatory instabilities are size dependent effects, some of them survive in the limit of an 
infinite lattice. Thus in the non-even potential case, both modes turn out to be unstable |34| . The recurrence relation 
(|13|l assures that SAB with frequencies slightly above the phonon band do not exist for B < (\/iKl/2 < I-K3Q. 
However, in this parameter region numerical simulations show that there exist breathers whose amplitudes do not 
tend to zero as wi, — > 2 + |60| | . As a consequence the energies of this family of breathesrs lay always above a certain 
positive lower bound. The existence of such an energy threshold is a rarely observed phenomenon in one-dimensional 
lattices. [43| 

The breathers obtained for B < have the same symmetries (site-centred and bond-centred modes) and the same 
stability properties as SAB with non-even potentials, i.e. both modes show oscillatory instabilities and also a harmonic 
instability in the case of the Sicvers-Takeno mode. 

Another interesting problem is the study of breathers superposed to uniformly stressed static states of the F PU chai n 
given by x n = cn, being c a nonzero constant |55| . Introducing the change of variables x n (t) = cn + x n (t \JW" '(c)), 
the previous analytical and numerical analysis can be completely applied to the renormalized displacements x n 
[60| . As a consequence, it is easy to prove that the FPU model has a family of exact breather solutions satisfying 
\mL n ^± 00 u n — c, which can be approximated by 



u n (t) ~ c+ (-W-J^ ^ , (19) 

W 1 ' V B(c) cosh(ny/ fi/W"(c)) V ' 

where n — u)\ — 4 W"(c) <C 1, and B(c) = |W"(c) (c) — (W^ (c)) 2 > 0. An example of a breather with a uniform 
stress is shown in Fig. with full circles. Note that the lower frequency of these breather solutions lies inside the 
phonon band if W"(c) < 1, and above it for W"(c) > 1. 



VII. DARK BREATHERS IN KLEIN-GORDON LATTICES 



The study of breathers in chains of nonlineraly interacting particles has a natural extension in the study of other 
localized entities called " dark breathers" . The term dark breather refers to a state of the chain where most of the 
oscillators are excited except one or a few units of them which have a very small amplitude. That is, in some sense 
it posses the opposite properties than a breather. This section summarizes some results about the conditions for the 
existence and stability properties of dark breathers in different Klein-Gordon lattices. 

Consider one-dimensional Klein-Gordon lattices of oscillators with Hamiltonian 

H = Y J ^l + v ( u n))+eW(u), (20) 

n 

where u n represents the coordinates of the oscillators referred to their equilibrium positions; V(u n ) represents the on- 
site potential; u represents the set of variables {u n }; and eW(u) represents the interacting potential. The parameter 
e describes the intensity of the coupling, and W{u) is given by 



W M= 2 EK+i-"") 2 - (21) 

n 

This interaction is attractive for e > 0, as a nonzero value of a variable tends to increase the values of the neighbouring 
variables with the same sign. The on-site potential is given by 

V(u n ) = ico 2 u 2 n + cf>(u n ), (22) 

with 4>(u n ) being the anharmonic part of the potential. The variables can be scaled so that all the particles in the 
lattice have mass unity and the linear frequency con = 1. 

The existence theorem by Mackay and Aubry [4(| establishes that dark breather solutions are possible, at least 
up to a certain value of the coupling parameter e c . Dark breathers can be calculated numerically. The method to 
calculate this type of solutions is similar to the used for obtaining breathers Fig. refFigl2, shows two different 
examples of dark breather profiles for a chain with a cubic soft on-site potential, with attractive interaction (left 
panel), and with repulsive interaction (right panel). 
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FIG. 12: Dark breathers profile for a cubic potential with attractive interaction (left) and with repulsive interaction(right), 
with e = 0.023. 



The linear stability analysis of dark breather solutions can be performed numerically. For that porpose both the 
Floquct multipliers and the Aubry band theory can be considered jlOl l50j . 

In the case of a soft on-site potential of the type (p(u n ) — ^\ u \ > an d attractive interaction between particles, dark 
breather solutions are not stable. Fig. ^|shows the type of stability of a dark breather and a breather by the Floquet 
multipliers method. The breather (left) is stable for e < 0.1, whereas the dark breather (right) is unstable for all the 
e values considered. This instability is due to harmonic bifurcations. 

Considering chains with repulsive interactions, the linear stability analysis shows that dark breather solutions are 
stable up to a significant value of the parameter e. 

The type of stability behaviour described above appears also in other kind of soft on-site potentials as, for example, 
the Morse potential given by 

V(u n ) = D (exp(-bu n ) - l) 2 . (23) 

Dark breather solutions are stable if the coupling is repulsive for e < 0.024. For bigger values the solution becomes 
unstable due to oscillatory and subharmonic bifurcations Fig. 1 141 shows the evolution of the Floquet multipliers 
for e = 0.024 (left panel), where the dark breather solution is still stable, and for e > 0.024 (right panel), where 
oscillatory and subharmonic bifurcations appear. 

Lattices with hard on-site potential present a different scenario. The dark breather solution is stable with attractive 
coupling potential, up to a certain value of e, and is unstable if the coupling is repulsive. Oscillatory, subharmonic 
and harmonic bifurcations appear depending on the form of the on-site potential. For a hard on-site potential of the 
type <p(u n ) = itt* and attractive coupling, the system is stable for e < 0.022. Figs. 1 151 and ^Jrepresent, respectively, 
the band structure and the Floquet multipliers at e = 0.041. They show that the dark breather solution is unstable 
due to an harmonic and small oscillatory bifurcations. The instability mode, shown in Fig. 1151 (right panel), is an 
asymmetric extended one. Simulations performed perturbing with it the dark breather give rise to a small oscillation 
with both sides of the chain out of phase, superimposed on the dark breather one, but the darkness is preserved. 
After that, there are oscillatory bifurcations due to the band mixing with the usual properties. 

An analysis of larger systems shows that subharmonic and oscillatory bifurcations persist even though the system 
is infinite, while in the case of breathers there are two kinds of size-dependent bifurcation |49| . 
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FIG. 13: Evolution of the Floquet multipliers with cubic on-site potential and attractive interaction for: (left) a breather with 
e — 0.1; (right) a dark breather with e = 0.004. The breather frequency is in both cases u>b = 0.8. 
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FIG. 14: Floquet eigenvalues for a Morse on-site potential and repulsive coupling. Left: e = 0.024, the system is still stable. 
Right: e = 0.035, the system becomes unstable due to oscillatory and subharmonic bifurcations. 

VIII. POLARONS IN BIOMOLECULES 

The study of charge transport (CT) in biomolecules, and particularly in DNA, is of utmost importance. Electronic 
transport plays a very important role in many biological processes, such as DNA repair after damage by radiation or 
chemical agents, photosynthesis, nevous conduction, etc. Also, because the quasi-one-dimensional structure of DNA 
and other biomolecules could be used to design electronic devices based on biomaterials (for a detailed bibliographical 
information, see i.e. [13 )• 

When an electron moves along the double strand of DNA it is accompanied by a local deformation of the molecule, 
forming a polaron or an electron-vibron breather. The polaronic character of the CT along DNA has been studied by 
means of different models, some of then based on the Peyrard-Bishop DNA model [53], in which the only variables 
are the distances between bases within each base pair, have permitted to show that the polaronic charge transport 
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FIG. 16: Left: Floquet multipliers with a quartic hard on-site potential and attractive interaction at e — 0.041 and ujb = 1.2. 
A harmonic and a small oscillatory instabilities appear. Right: the velocities components (the position ones are zero) of the 
eigenvalue corresponding to the harmonic instability. 



along DNA can be mediated by the coupling of the charge carrying unit with the hydrogen-bond deformations, 

In particular, it has been shown in some simple, but powerful models based on the proposed by Peyrard-Bishop 
[53| , that the polaronic CRT along DNA can be mediated by coupling of the charge carrying unit with the hydrogen- 
bond deformations, and that this transport survives to a small amount of parametric and structural disorder |37| . or 
by means of the coupling of this charge unit with the twist deformation of neighbor base pairs [S^ • In this last case, 
the moving polaron is slower than in the previous one, but it is more robust under the introduction of parametric 
disorder. These differents regimes apperars for different values of a parameter a, the couples the transference integral 
with the deformations of the hydrogen bonds, a parameter which at present there is not any realible experimental 
data. 

Consider a three-dimensional semi-classical tight -binding model for synthetically produced DNA polymers built up 
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FIG. 17: Sketch of the structure of the DNA model. The bases are represented by bullets. 



from single type of base pairs, i.e. either poly(dG)-poly(dC) or poly(dA)-poly(dT) DNA polymers Utilizing a 

nonlinear approach based on the concept of breather and polaron solutions it is possible to explore if the conductivity 
depends on the type of the DNA polymer. If they have quantitatively distinct transpor properties or not might be 
of interest for the design of synthetic molecular wires. Moreover, it has been possible to ameliorate preceding studies 
of electronic transfer in DNA in the sense that, instead of adjusting the coupling parameters, plausible 

estimates for them, derived with the help of quantum-chemical methods, has been used. 

The model for CT in DNA is based on the finding that the charge migration process is dominantly influenced by 
the transverse vibrations of the bases relative to each other in radial direction within a base pair plane. The impact 
of other vibrational degrees of freedom are expected to be negligible with respect to charge transport in DNA. Then, 
the motion can be viewed as confined to the base pair planes |6lj |. The Hamiltonian for CT in along a strand in 
DNA has two parts H = H e i + H V ib, with H e i describing the CT over the base pairs and H v n, the dynamics of radial 
vibrations of the base pairs or the vibronic part. A sketch of the structure of the DNA model is shown in Fig. 1171 

Supposing that vibrational dynamics can be described classically, it is easy to formulate the dynamical equations of 
the system (a detailed description of procedure can be found in 52]). It is possible to use some typical parameter values 
for DNA molecules and compute plausible values for the electron-coupling strengths as a result of quantum-chemical 
computational procedure [38|. 

The nonlinear interplay between the electronic and the vibrational degrees of freedom, makes possible the formation 
of polaronic electron-vibration compounds. Also, there exist some stationary solutions that can be calculated with 
the help of the nonlinear map approach explained in detail in |4^ . In general, the polarons are of fairly large 
extension (width). Regardless of the DNA polymer type, the electronic wave function is localized at a lattice site and 
the envelope of the amplitudes decays monotonically and exponentially as distance grows from this central site (base 
pair). However, the electronic wave function of the poly(dG)-poly(dC) DNA polymer is stronger localized than the 
one of its poly(dA)-poly(dT) counterpart. As polarons solutions are fairly large extended (the half-width involves 
< 100 lattice sites) it is plausible to think that the can be mobile. This polaron motion can be activated using the 
discrete gradient method to obtain suitable initial perturbations of the momentum coordinates which initiate coherent 
motion of the polaron compound. The electronic occupation probability is defined as n(i) = ^ n n |c„(i)| 2 , where the 
index n denotes the site of the n— th base on a strand and c„ determines the probability to find the electron (charge) 
residing at this site. 

Results relative to the time evolution of the first momentum of the electronic occupation probability are shown 
in Fig. ED Conductivity in synthetically produced DNA molecules depends on the type of the single base pair 
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of which the polymer is built of. While a polaron-like mechanism, relying on the nonlinear coupling between the 
electron amplitude and radial vibrations of the base pairs, is responsible for long-range and stable electronic transfer 
in (dG)-(dC) DNA polymers, the conductivity is comparatively weaker in the case of (dA)-(dT) DNA polymers. 
Especially when it comes to designing synthetic molecular wires these findings might be of interest. In fact, recent 
experiments suggest that electronic transfer through DNA molecules proceeds by polaron hopping |44| . Furthermore, 
these results comply with the findings of experiments which show also that poly(dG)-poly(dC) DNA polymers forms 
a better conductor than their poly(dA)-poly(dT) counterparts. 
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FIG. 18: Time evolution of the first momentum of the electronic occupation probability. Full (dashed) line: poly(dA)-poly(dT) 
(poly(dG)-poly(dC)) DNA polymer. 



IX. QUANTUM BREATHERS 

At present, the phenomenon of localization of energy due to nonlincarity in discrete classical lattices is relatively 
well understood, but the understanding of the quantum equivalence of discrete breathers is very limited. There exist 
some theoretical results (2^1 E3| and some experimental observations of these states in different quantum systems, 
as mixtures of 4-methyl-pyridine [3^|, in Cu benzoate 0, and in doped alkalihalides Also, the study of 

discrete breathers in quantum lattices supporting a small number of quanta is not only interesting theoretically, this 
knowledge may be relevant to studies of quantum dots and quantum computing |45| , and for studies of Bose-Einstcin 
condensates trapped in optical traps |2(. 

This section focus the attention in one-dimensional lattices with a small number of quanta described by the the 
quantum version of the discrete nonlinear Schrodinger equation (QDNLS), also known as the Hubbard model. This 
is a particularly simple model for a lattice of coupled anharmonic oscillators, which has been used to describe the 
dynamics of a great variety of systems [5]| ■ We present the results corresponding to a periodic lattice with / sites 
containing bosons |30j . Many of these results can be extended to a great variety of systems, for example, similar 
results have been obtained with a periodic lattice containing fermions, described by an attractive fermionic Hubbard 
model (FH) with two kinds of particles with opposite spins. This is a model of interest in connection with the theory 
of high-T c superconductivity, and it can be used to describe bound states of electron and holes in some nanostructures 
as nanorings (excitons) |29l| . 

The QDNLS Hamiltonian of the quantum system is 

/ i 

6 = - E 2^ b A bjbj + e $ ( ''< 1 + b ^ ^ 

3=1 
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where &j and bj are standard bosonic operators, Jj/cj is the ratio of anharmonicity to nearest neighbor hopping 
energy, and periodic boundary conditions are imposed. As first step, only short range interactions are considered 
given by the hopping coefficients ej. 

This Hamiltonian conmutes with the number operator N — X/f=i ^jbj, and can be block-diagonalised easily when 
the number of bosons is small enough. Numerically exact solutions can be found restricting the study to small lattices 
with a small number of quanta. The simplest nontrivial case is for N = 2, where bound states corresponds to bound 
two-vibron states, as has been observed experimentally in several systems |39j . These results have been extended to 
more complicated situations, noting that many of them are valid for larger values of N. 



A. Quantum breathers in a translational invariant lattice 



In a homogeneous quantum lattice the parameters jj and 6j are independent of j, i.e. jj — 7 and ej = e. With 
periodic boundary conditions, it is possible to block-diagonalize the Hamiltonian operator using eigenfunctions of the 
translation operator T, defined as Tfrj = bj +1 T. In each block, the eigenfunctions have a fixed value of the momentum 
k, with r = exp(ifc) being an eigenvalue of the translation operator |57j . In this way, it is possible to calculate the 
dispersion relation E(k) with a minimal computational effort. 

In general, as shown in ED if parameter 7 is high enough, there exists an isolated ground state eigenvalue for each k 
which corresponds to a quantum breather. In this state, there is a high probability of finding the quanta on the same 
site, but due to the translational invariance of the system, there exists an equal probability of finding these quanta 
at any site of the system. In these cases, some analytical expressions can be obtained in some asymptotic limits |57j . 




FIG. 19: Eigenvalues E{k). N = 2, / = 19 and 7 = 4. 



For example, using the number-state-basis \ip n >= [71°, nf, ...,n^], where nf represents the number of quanta at 
site i, in the case N = 2 and k = the (unnormalized) ground state is given by 

|*) = [20 ... 0] + [020 . . . 0] + • • • + [0 . . . 02] + 0(7 _1 ). 

The coefficients a, of the first / terms are equal to unity and the rest are 0(7 _1 ). 
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B. Trapping in lattices with broken translational symmetry 



When the translational invariance of a lattice is broken, the Hamiltonian operator does not commute with the 
translation operator and breathers can be localized around a particular point of the lattice. Perhaps the simplest way 
to do it is to consider non-flux boundary conditions, in order to simulate a finete-size chain. In this case, the ground 
state becomes localized around the middle of the chain but, when / is high enough, this effect vanishes. 

The existence of local inhomogeneities or impurities can also break the translational symmetry of the system, and 
this can be modeled by making some coefficients dependent on the site. This can occurs in system with twisted or 
bent geometries, as in models of globular proteins [3l| or photonic crystals In these cases, the geometry effects can 
be considered introducing a long-range hopping term in the Hamiltonian operator. 

The breaking of translational symmetry of the system localizes the breather state around a particular site of the 
chain, as shown in Fig. |23 Note that the localization phenomena due to random variation of lattice parameters have 
been studied in harmonic models since the work of Anderson Q but, in all cases, and due to anharmonic effects, 
there exist new localization effects, and, when it is present, the anharmonicity enhances the Anderson-like effect. 
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FIG. 20: Square wave function amplitudes corresponding to localized components of the ground state. / = 11, TV = 4 and 
7 = 2. Point impurity at the site I = 6. (• • •) homogeneous chain case, 7i m = 2, ( ) 7i m = 2.1 , ( — ) 7i m = 2.5. 



X. THERMAL EVOLUTION OF ENZYME-CREATED OSCILLATING BUBBLES IN DNA 

This last section reviews some fundamental researches on the creation and subsequent evolution of breathers in DNA 
chains. For a double-stranded DNA molecule to replicate, the two strands of the double helix must be separated from 
each other, at least locally, with the breaking up of the hydrogen bonds. For long double-stranded DNA molecules, 
the rate of spontaneous strand separation is very low under physiological condition. Specific enzimes belonging to the 
large family of nucleic acid helicases (see jg^) (there are more than sixty different types), utilize the energy of ATP 
hydrolysis to power strand separation. The opening of a selected region of the twisted (Watson-Crick) DNA double 
helix, begins with a partial unwinding at an area called the replication fork and is observed as a bubble. Helicases 
are ubiquitous and versatile enzymes that, in conjunction with other components of the macromolecular machines, 
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carry out important biological processes such as DNA replication; DNA recombination; nucleotide excision repair 
(i.e. UV damage); RNA editing and splicing; the transfer of a single-stranded nucleic acid to other nucleic acid or a 
protein or its release into solution Q2J and formation of T-loop telomeres 63]. Nucleic acid helicases have several 
structural varieties (monomeric, dimeric, trimeric, tetrameric and closed hcxamcric), but all of them use the hydrolysis 
of nucleoside triphosphate (NTP) to nucleoside diphosphate (NDP) as the preferred source of energy. Helicases have 
also a huge interest from a point of view of clinical disease. Actually, there are known up to five human monogenic 
hereditary disorders in which the mutated protein responsible for the disease carries an aberrant helicase activity: 
cockayne 's syndrome, xeroderma pigmentosum, trichothiodystrophy, Bloom 's syndrome, and Werner 's syndrome. 
The study of the formation of DNA openings has been initiated recently utilizing models based on nonlinear lattice 
dynamics Q, . 

With a view to vibrational excitations in DNA, the Peyrard-Bishop (PB) model [5^| and its successors [Til l2r| 
have been successfully applied to describe moving localized excitations (moving breathers), which reproduce typical 
features of the DNA opening dynamics such as the magnitude of the amplitudes and the time scale of the " oscillating 
bubble" preceding full strand separation. 

The bubble formation process is initiated by structural deformations of selected regions of the parental DNA duplex 
that serves as a template for the replication. The process begins when the replication apparatus identifies the starting 
point and then gets combined to it. The replication apparatus is and assembled complex that forms the replication fork 
and opens it directionally. During this process the helical DNA is unwound and replicated. One of the components of 
that complex a type of helicase, which is unable to recognize the origin of the replication per se, and which requires 
the participation of specific proteins to lead it to the initiation site |27|. 

It is assumed that this enzyme operates in the way, that a local unwinding in combination with (rather small) 
stretchings of the H-bonds in this region occurs. A working model should demonstrate that there are indeed initial 
deformations that give rise to localized vibrations, constrained to a region of DNA, matching the properties of oscil- 
lating bubbles observed experimentally. These oscillating bubbles with their temporarily extended but yet unbroken 
H-bonds serve as the precursors to the replication bubble. Furthermore, it is important to know whether the stable 
radial and torsional breathers persist under imposed thermal perturbations. 

The research has been done using a nonlinear oscillator network model for the DNA double helix that is explained in 
detail in 0, IT3. ITsj . The equilibrium position of each base within the duplex configuration is described in a Cartesian 

coordinate system by xf\, \ and z^\. The index pair (n,i) labels the n-th base on the i— th strand with i = 1,2 
and 1 ^ n ^ N , where N is the number of base pairs considered. Displacements of the bases from their equilibrium 
positions are denoted by x n ,i, Dn,i an d 2n,i- The potential energy taking into account the interactions between the 
bases consists of four parts. The potential energy of the hydrogen bond within a base pair is modeled typically by a 
Morse potential 



v n — n 



exp 



(-f d n ) -l] , (25) 



where the variables d n describe dynamical deviations of the hydrogen bonds from their equilibrium lengths do- The 
site-dependent depth of the Morse potential, D n , depends on the number of involved hydrogen bonds for the two 
different pairings in DNA, namely the G-C and the A-T pairs. The former pair includes three hydrogen bonds while 
the latter includes only two. a -1 is a measure of the potential-well width. 

The energies of the rather strong and rigid covalent bonds between the nucleotides n and n — 1 on the i — th strand 
are modeled by harmonic potential terms 

vr = Y l l*> ( 26 ) 

and l n ,i describes the deviations from the equilibrium distance between two adjacent bases on the same strand. K is 
the elasticity coefficient. 

Effects of stacking, which impede that, due to the backbone rigidity, one base slides over another [(5l| are incorpo- 
rated in the following potential terms 

Vr = - (d n>i - d„_ M ) 2 . (27) 

The supposedly small deformations in longitudinal direction can be modeled by harmonic elasticity potential terms 
given by 



Vf'^ji^i-Zn-uf. (2M 
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FIG. 21: The initial radial distortions of the double helix. The inter-base distance d n (0) in 



The kinetic energy of a nucleotide is determined by 



El 



1 

2m 



OS) 2 + (pS) 2 + (pI: 



(29) 



where m is the mass and p^ f ^ denotes the (x,y, z)— component of the momentum. 



The model Hamiltonian reads then as 



JY 



H 



EE E " 

t=i,2 n=l 



with 



(30) 



(31) 



and the summation in H30(l is performed over all nucleotides and the two strands. 

Consider the initiation of the oscillating bubbles in DNA. The starting point is a DNA molecule for which a certain 
segment experiences initially angular and radial deformations due to the action of some helicase to which a region 
of the DNA is bound. In order to simulate the deforming action of enzymes, assume that initially a number of 
consecutive sites in the center of the DNA lattice (hereafter referred to as the central region) are exerted to forces 
acting in angular and radial direction such that in this region the molecule experiences twist reduction together with 
radial stretchings. These structural deformations can be extended over a region encompassing up to thirty base pairs 
and as can be demonstrated give rise to the formation of H-bridge breather solutions (ex tending over 15—20 base 
pairs) reproducing the "oscillating bubbles" observed for the DNA-opening process P3 • In Figs. El and 01 the 
localized initial distortions are shown. Both the angular and radial deformation patterns are bell-shaped, due to the 
fact that the enzymatic force is exerted locally, i.e. in the extreme case to a single base pair only for which the H-bond 
is deformed 

The first one being of non-positive amplitudes is linked with reduced twist while the latter one with non-negative 
amplitudes is associated with radial stretchings. The radial and angular deformation patterns are centered at the 
central lattice site (base pair) at which the H-bond stretching and twist angle reduction is at maximum. At either 
side of the central site the amplitudes approach progressively zero. The deformation energy amounts to 0.0362 eV. 
The associated set of coupled equations has been integrated with the help of a fourth-order Runge-Kutta method. 
For the simulation the DNA lattice consists of 799 sites and open boundary conditions are imposed. The same initial 
conditions are used both at zero and nonzero temperatures. 

The study of the temperature dependence of these oscillating bubbles has been done by means of bringing in contact 
the nonlinear oscillator network with a heat bath using the Nos-Hoover-method. It is demonstrated that the radial and 
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FIG. 22: The initial angular distortions of the double helix. 

torsional breathers sustain the impact of thermal perturbations, even at temperatures as high as room temperature. 
For T^O the breathers start to travel coherently towards an end of te DNA duplex, representing realistically the 
oscillating bubbles found in NAs i.e. they oscillate with periods in the range 0.3 -0.8 ps, are on the average extended 
over 10 -20 base pairs and possess maximal amplitudes of the order of 0.6 A. 

Quantitatively, the results regarding the temperature dependence of the breather evolution are suitably summarized 
by the time-evolution of the first momentum of the energy distribution defined as 

N 

n(t) =J2 ^{n c -n) E n ^{t) , (32) 

1=1,2 n=l 

and the energy E n .i is defined in Eq. I|31|l and n c is the site index corresponding to the center of the DNA lattice. 

This quantity describes the temporal behavior of the position of the center of a breather. Thus it represents 
a measure for the mobility of the breathers. Generally, the higher the temperature the larger the amplitude of the 
radial breather becomes and the faster the radial (and torsional) breather travels along the DNA lattice. This behavior 
is illustrated in Fig. [23 

In summary, out of an initial non-equilibrium situation, for which the hydrogen bonds of an under-twisted segment 
of the DNA lattice have been stretched, breathers develop in the radial and angular displacement variables. For zero 
temperature the breathers remain standing in the initially excited region. However, for T > the breathers start to 
travel coherently towards an end of the DNA duplex. The observed breathers represent realistically the oscillating 
bubbles found prior to complete unzipping of DNA, i.e. they oscillate with periods in the range 0.3 — O.Sps, are on 
the average extended over 10 — 20 base pairs and possess maximal amplitudes of the order of < 0.6 A. These results 
demonstrate that the action of some enzyme, mimicked by localized radial and torsional distortions of the DNA 
equilibrium configuration, initiates in fact the production of oscillating bubbles in DNA. Moreover, these "oscillating 
bubbles" sustain the impact of thermal perturbations. 
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